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Abstract 

Cooperative diversity has been recently proposed as a way to form virtual antenna arrays that provide dramatic 
gains in slow fading wireless environments. However most of the proposed solutions require distributed space-time 
coding algorithms, the careful design of which is left for future investigation if there is more than one cooperative 
relay. We propose a novel scheme, that alleviates these problems and provides diversity gains on the order of the 
number of relays in the network. Our scheme first selects the best relay from a set of M available relays and 
then uses this "best" relay for cooperation between the source and the destination. We develop and analyze a 
distributed method to select the best relay that requires no topology information and is based on local measurements 
of the instantaneous channel conditions. This method also requires no explicit communication among the relays. 
The success (or failure) to select the best available path depends on the statistics of the wireless channel, and a 
methodology to evaluate performance for any kind of wireless channel statistics, is provided. Information theoretic 
analysis of outage probability shows that our scheme achieves the same diversity-multiplexing tradeoff as achieved 
by more complex protocols, where coordination and distributed space-time coding for M nodes is required, such 
as those proposed in [7]. The simplicity of the technique, allows for immediate implementation in existing radio 
hardware and its adoption could provide for improved flexibility, reliability and efficiency in future 4G wireless 
systems. 

Index Terms 

Network cooperative diversity, outage probability, coherence time, fading channel, wireless networks. 

I. Introduction 

IN this work, we propose and analyze a practical scheme that forms a virtual antenna array among single 
antenna terminals, distributed in space. The setup includes a set of cooperating relays which are willing to 
forward received information towards the destination and the proposed method is about a distributed algorithm that 
selects the most appropriate relay to forward information towards the receiver. The decision is based on the end- 
to-end instantaneous wireless channel conditions and the algorithm is distributed among the cooperating wireless 
terminals. 

The best relay selection algorithm lends itself naturally into cooperative diversity protocols [12], [13], [6], [14], 
which have been recently proposed to improve reliability in wireless communication systems using distributed virtual 
antennas. The key idea behind these protocols is to create additional paths between the source and destination using 
intermediate relay nodes. In particular, Sendonaris, Erkip and Aazhang [12], [13] proposed a way of beamforming 
where source and a cooperating relay, assuming knowledge of the forward channel, adjust the phase of their 
transmissions so that the two copies can add coherently at the destination. Beamforming requires considerable 
modifications to existing RF front ends that increase complexity and cost. Laneman, Tse and Wornell [6] assumed no 
CSI at the transmitters and therefore assumed no beamforming capabilities and proposed the analysis of cooperative 
diversity protocols under the framework of diversity-multiplexing tradeoffs. Their basic setup included one sender, 
one receiver and one intermediate relay node and both analog as well as digital processing at the relay node 
were considered. Subsequently, the diversity-multiplexing tradeoff of cooperative diversity protocols with multiple 
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Opportunistic Relaying 

Fig. 1. A transmission is overheard by neighboring nodes. Distributed Space-Time coding is needed so that all overhearing nodes could 
simultaneously transmit. In this work we analyze "Opportunistic Relaying" where the relay with the strongest transmitter-relay-receiver path 
is selected, among several candidates, in a distributed fashion using instantaneous channel measurements. 



relays was studied in [7], [1]. While [7] considered the case of orthogonal transmission between source and 
relays, [1] considered the case where source and relays could transmit simultaneously. It was shown in [1] that 
by relaxing the orthogonality constraint, a considerable improvement in performance could be achieved, albeit at a 
higher complexity at the decoder. These approaches were however information theoretic in nature and the design 
of practical codes that approach these limits was left for further investigation. 

Such code design is difficult in practice and an open area of research: while space time codes for the Multiple 
Input Multiple Output (MIMO) link do exist [16] (where the antennas belong to the same central terminal), more 
work is needed to use such algorithms in the relay channel, where antennas belong to different terminals distributed 
in space. The relay channel is fundamentally different than the point-to-point MIMO link since information is 
not a priori known to the cooperating relays but rather needs to be communicated over noisy links. Moreover, 
the number of participating antennas is not fixed since it depends on how many relay terminals participate and 
how many of them are indeed useful in relaying the information transmitted from the source. For example, for 
relays that decode and forward, it is necessary to decode successfully before retransmitting. For relays that amplify 
and forward, it is important to have a good received SNR, otherwise they would forward mostly their own noise 
[20]. Therefore, the number of participating antennas in cooperative diversity schemes is in general random and 
space-time coding invented for fixed number of antennas should be appropriately modified. It can be argued that 
for the case of orthogonal transmission studied in the present work (i.e. transmission during orthogonal time or 
frequency channels) codes can be found that maintain orthogonality in the absence of a number of antennas (relays). 
That was pointed in [7] where it was also emphasized that it remains to be seen how such codes could provide 
residual diversity without sacrifice of the achievable rates. In short, providing for practical space-time codes for 
the cooperative relay channel is fundamentally different than space-time coding for the MIMO link channel and is 
still an open and challenging area of research. 

Apart from practical space-time coding for the cooperative relay channel, the formation of virtual antenna 
arrays using individual terminals distributed in space, requires significant amount of coordination. Specifically, 
the formation of cooperating groups of terminals involves distributed algorithms [7] while synchronization at the 
packet level is required among several different transmitters. Those additional requirements for cooperative diversity 
demand significant modifications to almost all layers of the communication stack (up to the routing layer) which 
has been built according to "traditional", point-to-point (non-cooperative) communication. 

In fig.n a transmitter transmits its information towards the receiver while all the neighboring nodes are in listening 
mode. For a practical cooperative diversity in a three-node setup, the transmitter should know that allowing a relay 
at location (B) to relay information, would be more efficient than repetition from the transmitter itself. This is not 

'Note that in that scheme the relays do not transmit in mutually orthogonal time/frequency bands. Instead they use a space-time code to 
collaboratively send the message to the destination. Orthogonality refers to the fact that the source transmits in time slots orthogonal to the 
relays. Throughout this paper we will refer to Laneman's scheme as orthogonal cooperative diversity. 
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a trivial task and such event depends on the wireless channel conditions between transmitter and receiver as well 
as between transmitter-relay and relay-receiver. What if the relay is located in position (A)? This problem also 
manifests in the multiple relay case, when one attempts to simplify the physical layer protocol by choosing the 
best available relay. In [19] it was suggested that the best relay be selected based on location information with 
respect to source and destination based on ideas from geographical routing proposed in [17]. Such schemes require 
knowledge or estimation of distances between all relays and destination and therefore require either a) infrastructure 
for distance estimation (for example GPS receivers at each terminal) or b) distance estimation using expected SNRs 
which is itself a non-trivial problem and is more appropriate for static networks and less appropriate for mobile 
networks, since in the latter case, estimation should be repeated with substantial overhead. 

In contrast, we propose a novel scheme that selects the best relay between source and destination based on 
instantaneous channel measurements. The proposed scheme requires no knowledge of the topology or its estimation. 
The technique is based on signal strength measurements rather than distance and requires a small fraction of the 
channel coherence time. All these features make the design of such a scheme highly challenging and the proposed 
solution becomes non-trivial. Additionally, the algorithm itself provides for the necessary coordination in time and 
group formation among the cooperating terminals. 

The three-node reduction of the multiple relay problem we consider, greatly simplifies the physical layer design. 
In particular, the requirement of space-time codes is completely eliminated if the source and relay transmit in 
orthogonal time-slots. We further show that there is essentially no loss in performance in terms of the diversity- 
multiplexing tradeoff as compared to the transmission scheme in [7] which requires space-time coding across the 
relays successful in decoding the source message. We also note that our scheme can be used to simplify the non- 
orthogonal multiple relay protocols studied in [1]. Intuitively, the gains in cooperative diversity do not come from 
using complex schemes, but rather from the fact that we have enough relays in the system to provide sufficient 
diversity. 

The simplicity of the technique, allows for immediate implementation in existing radio hardware. An implementa- 
tion of the scheme using custom radio hardware is reported in [3]. Its adoption could provide for improved flexibility 
(since the technique addresses coordination issues), reliability and efficiency (since the technique inherently builds 
upon diversity) in future 4G wireless systems. 

A. Key Contributions 

One of the key contribution of this paper is to propose and analyze a simplification of user cooperation protocols 
at the physical layer by using a smart relay selection algorithm at the network layer. Towards this end, we take the 
following steps: 

• A new protocol for selection of the "best" relay between the source and destination is suggested and analyzed. 
This protocol has the following features: 

- The protocol is distributed and each relay only makes local channel measurements. 

- Relay selection is based on instantaneous channel conditions in slow fading wireless environments. No 
prior knowledge of topology or estimation of it is required. 

- The amount of overhead involved in selecting the best relay is minimal. It is shown that there is a flexible 
tradeoff between the time incurred in the protocol and the resulting error probability. 

• The impact of smart relaying on the performance of user cooperation protocols is studied. In particular, it 
is shown that for orthogonal cooperative diversity protocols there is no loss in performance (in terms of the 
diversity multiplexing tradeoff) if only the best relay participates in cooperation. Opportunistic relaying provides 
an alternative solution with a very simple physical layer to conventional cooperative diversity protocols that 
rely on space-time codes. The scheme could be further used to simplify space-time coding in the case of 
non-orthogonal transmissions. 

Since the communication scheme exploits the wireless channel at its best, via distributed cooperating relays, 
we naturally called it opportunistic relaying. The term "opportunistic" has been widely used in various different 
contexts. In [23], it was used in the context of repetitive transmission of the same information over several paths, 
in 802. 1 lb networks. In our setup, we do not allow repetition since we are interested in providing diversity without 
sacrificing the achievable rates, which is a characteristic of repetition schemes. The term "opportunistic" has also 
been used in the context of efficient flooding of signals in multi-hop networks [24], to increase communication 
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Fig. 2. Source transmits to destination and neighboring nodes overhear the communication. The "best" relay among M candidates is selected 
to relay information, via a distributed mechanism and based on instantaneous end-to-end channel conditions. For the diversity-multiplexing 
tradeoff analysis, transmission of source and "best" relay occur in orthogonal time channels. The scheme could be easily modified to 
incorporate simultaneous transmissions from source and "best" relay. 



range and therefore has no relationship with our work. We first encountered the term "opportunistic" in the work 
by Viswanath, Tse and Laroia [25], where the base station always selects the best user for transmission in an 
artificially induced fast fading environment. In our work, a mechanism of multi-user diversity is provided for 
the relay channel, in single antenna terminals. Our proposed scheme, resembles selection diversity that has been 
proposed for centralized multi-antenna receivers [8]. In our setup, the single antenna relays are distributed in space 
and attention has been given in selecting the "best" possible antenna, well before the channel changes again, using 
minimal communication overhead. 

In the following section, we describe in detail opportunistic relaying and contrast its distributed, location information- 
free nature to existing approaches in the field. Probabilistic analysis and close form expressions regarding the success 
(or failure) and speed of "best" path selection, for any kind of wireless channel statistics, are provided in section 
ITTTl In section |W] we prove that opportunistic relaying has no performance loss compared to complex space-time 
coding, under the same assumptions of orthogonal channel transmissions [7] and discuss the ability of the scheme 
to further simplify space-time coding for non-orthogonal channels. We also discuss in more detail, why space-time 
codes designed for the MIMO link are not directly applicable to the cooperative relay channel. We conclude in 
section [V] 

II. Description of Opportunistic Relaying 

According to opportunistic relaying, a single relay among a set of M relay nodes is selected, depending on which 
relay provides for the "best" end-to-end path between source and destination (fig. ^ GJ- The wireless channel a s i 
between source and each relay i, as well as the channel between relay i and destination affect performance. These 
parameters model the propagation environment between any communicating terminals and change over time, with 
a rate that macroscopically can be modelled as the Doppler shift, inversely proportional to the channel coherence 
time. Opportunistic selection of the "best" available relay involves the discovery of the most appropriate relay, in 
a distributed and "quick" fashion, well before the channel changes again. We will explicitly quantify the speed of 
relay selection in the following section. 

The important point to make here is that under the proposed scheme, the relay nodes monitor the instantaneous 
channel conditions towards source and destination, and decide in a distributed fashion which one has the strongest 
path for information relaying, well before the channel changes again. In that way, topology information at the relays 
(specifically location coordinates of source and destination at each relay) is not needed. The selection process reacts 
to the physics of wireless propagation, which are in general dependent on several parameters including mobility 
and distance. By having the network select the relay with the strongest end-to-end path, macroscopic features 
like "distance" are also taken into account. Moreover, the proposed technique is advantageous over techniques 
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that select the best relay a priori, based on distance toward source or destination, since distance-dependent relay 
selection neglects well-understood phenomena in wireless propagation such as shadowing or fading: communicating 
transmitter-receiver pairs with similar distances might have enormous differences in terms of received SNRs. 
Furthermore, average channel conditions might be less appropriate for mobile terminals than static. Selecting 
the best available path under such conditions (zero topology information, "fast" relay selection well bellow the 
coherence time of the channel and minimum communication overhead) becomes non-obvious and it is one of the 
main contributions of this work. 

More specifically, the relays overhear a single transmission of a Ready-to-Send (RTS) packet and a Clear-to- 
Send (CTS) packet from the destination. From these packets the relays assess how appropriate each of them is 
for information relaying. The transmission of RTS from the source allows for the estimation of the instantaneous 
wireless channel a s i between source and relay i, at each relay i (fig. Similarly, the transmission of CTS from 
the destination, allows for the estimation of the instantaneous wireless channel between relay i and destination, 
at each relay i, according to the reciprocity theorem[26] 2 . Note that the source does not need to listen to the CTS 
packet 3 from the destination. 

Since communication among all relays should be minimized for reduced overall overhead, a method based on 
time was selected: as soon as each relay receives the CTS packet, it starts a timer from a parameter hi based on the 
instantaneous channel measurements a s i,aid- The timer of the relay with the best end-to-end channel conditions 
will expire first. That relay transmits a short duration flag packet, signaling its presence. All relays, while waiting 
for their timer to reduce to zero (i.e. to expire) are in listening mode. As soon as they hear another relay to flag 
its presence or forward information (the best relay), they back off. 

For the case where all relays can listen source and destination, but they are "hidden" from each other (i.e. they 
can not listen each other), the best relay notifies the destination with a short duration flag packet and the destination 
notifies all relays with a short broadcast message. 

The channel estimates a s i, at each relay, describe the quality of the wireless path between source-relay- 
destination, for each relay i. Since the two hops are both important for end-to-end performance, each relay should 
quantify its appropriateness as an active relay, using a function that involves the link quality of both hops. Two 
functions are used in this work: under policy I, the minimum of the two is selected (equation Q), while under 
policy II, the harmonic mean of the two is used (equation ©). Policy I selects the "bottleneck" of the two paths 
while Policy II balances the two link strengths and it is a smoother version of the first one. 

Under policy I: , • n 12 1 121 , n 

hi = mm{|a si | , \a id \ } (1) 



Under policy II: i 2 2 \a s i\ 2 la^' 2 

+ ~~ Kil 2 + 



hi — 1 , 1 — I 12 1 1- @) 



The relay i that maximizes function hi is the one with the "best" end-to-end path between initial source and final 
destination. After receiving the CTS packet, each relay i will start its own timer with an initial value T%, inversely 
proportional to the end-to-end channel quality hi, accjprding\ to the following equation: ^ 

K 

Here A is a constant. The units of A depend on the units of hi. Since hi is a scalar, A has the units of time. For 
the discussion in this work, A has simply values of usees. 

h b = max{/ij}, <S=^- (4) 
T b = min{Tj, i £ [1..M]. (5) 

Therefore, the "best" relay has its timer reduced to zero first (since it started from a smaller initial value, according 
to equations I©-©. This is the relay b that participates in forwarding information from the source. The rest of the 
relays, will overhear the "flag" packet from the best relay (or the destination, in the case of hidden relays) and 
back off. 

2 We assume that the forward and backward channels between the relay and destination are the same from the reciprocity theorem. Note 
that these transmissions occur on the same frequency band and same coherence interval. 

3 The CTS packet name is motivated by existing MAC protocols. However unlike the existing MAC protocols,the source does not need to 
receive this packet. 
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After the best relay has been selected, then it can be used to forward information towards the destination. 
Whether that "best" relay will transmit simultaneously with the source or not, is completely irrelevant to the relay 
selection process. However, in the diversity-multiplexing tradeoff analysis in section Hvl we strictly allow only one 
transmission at each time and therefore we can view the overall scheme as a two-step transmission: one from source 
and one from "best" relay, during a subsequent (orthogonal) time channel (fig. 

A. A note on Time Synchronization 

In principle, the RTS/CTS transmissions between source and destination, existent in many Medium Access Control 
(MAC) protocols, is only needed so that all intermediate relays can assess their connectivity paths towards source 
and destination. The reception of the CTS packet triggers at each relay the initiation of the timing process, within 
an uncertainty interval that depends on different propagation times, identified in detail in section Therefore, an 
explicit time synchronization protocol among the relays is not required. Explicit time synchronization would be 
needed between source and destination, only if there was no direct link between them. In that case, the destination 
could not respond with a CTS to a RTS packet from the source, and therefore source and destination would need 
to schedule their RTS/CTS exchange by other means. In such cases "crude" time synchronization would be useful. 
Accurate synchronization schemes, centralized [2] or decentralized [4], do exist and have been studied elsewhere. 
We will assume that source and destination are in communication range and therefore no synchronization protocols 
are needed. 

B. A note on Channel State Information ( CSI) 

CSI at the relays, in the form of link strengths (not signal phases), is used at the network layer for "best" relay 
selection. CSI is not required at the physical layer and is not exploited either at the source or the relays. The 
wireless terminals in this work do not exploit CSI for beamforming and do not adapt their transmission rate to 
the wireless channel conditions, either because they are operating in the minimum possible rate or because their 
hardware do not allow multiple rates. We will emphasize again that no CSI at the physical layer is exploited at the 
source or the relays, during the diversity-multiplexing tradeoff analysis, in section Hvl 

C. Comparison with geometric approaches 

As can be seen from the above equations, the scheme depends on the instantaneous channel realizations or 
equivalently, on received instantaneous SNRs, at each relay. An alternative approach would be to have the source 
know the location of the destination and propagate that information, alongside with its own location information to 
the relays, using a simple packet that contained that location information. Then, each relay, assuming knowledge of 
its own location information, could assess its proximity towards source and destination and based on that proximity, 
contend for the channel with the rest of the relays. That is an idea, proposed by Zorzi and Rao [17] in the context 
of fading-free wireless networks, when nodes know their location and the location of their destination (for example 
they are equipped with GPS receivers). The objective there was to study geographical routing and study the average 
number of hops needed under such schemes. All relays are partitioned into a specific number of geographical regions 
between source and destination and each relay identifies its region using knowledge of its location and the location 
of source and destination. Relays at the region closer to the destination contend for the channel first using a standard 
CSMA splitting scheme. If no relays are found, then relays at the second closest region contend and so on, until 
all regions are covered, with a typical number of regions close to 4. The latency of the above distance-dependent 
contention resolution scheme was analyzed in [18]. 

Zorzi and Rao's scheme of distance-dependent relay selection was employed in the context of Hybrid- ARQ, 
proposed by Zhao and Valenti [19]. In that work, the request to an Automatic Repeat Request (ARQ) is served 
by the relay closest to the destination, among those that have decoded the message. In that case, code combining 
is assumed that exploits the direct and relayed transmission (that's why the term Hybrid was used) 4 . Relays are 
assumed to know their distances to the destination (valid for GPS equipped terminals) or estimate their distances 
by measuring the expected channel conditions using the ARQ requests from the destination or using other means. 

4 The idea of having a relay terminal respond to an ARQ instead of the original source, was also reported and analyzed in [6] albeit for 
repetition coding instead of hybrid code combining. 
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Fig. 3. The middle row corresponds to the "best" relay. Other relays (top or bottom row) could erroneously be selected as "best" relays, 
if their timer expired within intervals when they can not hear the best relay transmission. That can happen in the interval [tL,tc] for case 
(a) (No Hidden Relays) or \tL,tu] for case (b) (Hidden Relays). t^,tj are time points where reception of the CTS packet is completed at 
best relay b and relay j respectively. 



We note that our scheme of opportunistic relaying differs from the above scheme in the following aspects: 

• The above scheme performs relay selection based on geographical regions while our scheme performs selection 
based on instantaneous channel conditions. In wireless environment, the latter choice could be more suitable 
as relay nodes located at similar distance to the destination could have vastly different channel gains due to 
effects such as fading. 

• The above scheme requires measurements to be only performed once if there is no mobility among nodes but 
requires several rounds of packet exchanges to determine the average SNR. On the other hand opportunistic 
relaying requires only three packet exchanges in total to determine the instantaneous SNR, but requires that 
these measurements be repeated in each coherence interval. We show in the subsequent section that the overhead 
of relay selection is a small fraction of the coherence interval with collision probability less than 0.6%. 

• We also note that our protocol is a proactive protocol since it selects the best relay before transmission. The 
protocol can easily be made to be reactive (similar to [19]) by selecting the relay after the first phase. However 
this modification would require all relays to listen to the source transmission which can be energy inefficient 
from a network sense. 

III. Probabilistic Analysis of Opportunistic Relaying 

The probability of having two or more relay timers expire "at the same time" is zero. However, the probability 
of having two or more relay timers expire within the same time interval c is non zero and can be analytically 
evaluated, given knowledge of the wireless channel statistics. 

The only case where opportunistic relay selection fails is when one relay can not detect that another relay is 
more appropriate for information forwarding. Note that we have already assumed that all relays can listen initial 
source and destination, otherwise they do not participate in the scheme. We will assume two extreme cases: a) all 
relays can listen to each other b) all relays are hidden from each other (but they can listen source and destination). 
In that case, the flag packet sent by the best relay is received from the destination which responds with a short 
broadcast packet to all relays. Alternatively, other schemes based on "busy tone" (secondary frequency) control 
channels could be used, requiring no broadcast packet from the destination and partly alleviating the "hidden" 
relays problem. 

From fig. |3j collision of two or more relays can happen if the best relay timer Tj, and one or more other relay 
timers expire within [tL, tc] for the case of no hidden relays (case (a)). This interval depends on the radio switch 
time from receive to transmit mode d s and the propagation times needed for signals to travel in the wireless 
medium. In custom low-cost transceiver hardware, this switch time is typically on the order of a few /^seconds 
while propagation times for a range of 100 meters is on the order of 1/3 fiseconds. For the case of "hidden" relays 
the uncertainty interval becomes [tL, tu] since now the duration of the flag packet should be taken into account, 
as well as the propagation time towards the destination and back towards the relays and the radio switch time at 
the destination. The duration of the flag packet can be made small, even one bit transmission could suffice. In any 
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case, the higher this uncertainty interval, the higher the probability of two or more relay timers to expire within that 
interval. That's why we will assume maximum values of c, so that we can assess worst case scenario performance. 

(a) No Hidden Relays: 

C — fmax ~\~ y^b f^j\max ~i~ d>s (6) 

(b) Hidden Relays: 

C ^max ~t~ 1^6 Tlj\max ~t~ -\- dlLT -\- 2jTl 7nax (7) 

• nf propagation delay between relay j and destination. n max is the maximum. 

• r: propagation delay between two relays. r max is the maximum. 

• d s : receive-to-transmit switch time of each radio. 

• dur: duration of flag packet, transmitted by "best" relay. 

In any case, the probability of having two or more relays expire within the same interval c, out of a collection of 
M relays, can be described by the following expression: 

Pr(Collision) < Pr(any Tj < T b + c | j / b) (8) 
where T b = min{Tj}, j S [1,M] and c > 0. 

Notice that we assume failure of relay selection when two or more relays collide. Traditional CSMA protocols 
would require the relays to sense that collision, back-off and retry. In that way collision probability could be further 
reduced, at the expense of increased latency overhead for relay selection. We will analyze the collision probability 
without any contention resolution protocol and further improvements are left for future work. 

We will provide an analytic way to calculate a close-form expression of equation © for any kind of wireless 
fading statistics. But before doing so, we can easily show that this probability can be made arbitrary small, close 
to zero. 

If T b = mm{Tj},j G [1, M] and Yi < Y 2 < . . . < Y M the ordered random variables {Tj} with T b = Y\, and Y 2 
the second minimum timer, thenp r((my Tj K Tb + c \ j ^ b) = Pr(Y 2 <Y 1 + c) (9) 

From the last equation, we can see that this probability can be made arbitrarily small by decreasing the parameter 
c. For short range radios (on the order of 100 meters), this is primarily equivalent to selecting radios with small 
switch times (from receive to transmit mode) on the order of a few microseconds. 

Given that Yj = X/hu), Y\ < Y 2 < . . . < Ym is equivalent to l/hm < l//&(2) < ■ • • < 1/^(m) 5 > equation © is 

equivalent t0 Pr(Y 2 < Y 1 + C ) = Pr(-i- < J- + V (10) 

h(2) h(i) A 

and Y± <Y 2 < ... < Ym O hm > hr 2 ) . . . > (h, A, c are positive numbers). 

From the last equation dlOL it is obvious that increasing A at each relay (in equation @), reduces the probability 
of collision to zero since equation dlOt goes to zero with increasing A. 

In practice, A can not be made arbitrarily large, since it also "regulates" the expected time, needed for the network 
to find out the "best" relay. From equation © and Jensen's inequality we can see that 

E[Tj) = E[\/hj] > X/E[hj] (11) 

or in other words, the expected time needed for each relay to flag its presence, is lower bounded by A times a 
constant. Therefore, there is a tradeoff between probability of collision and speed of relay selection. We need to 
have A as big as possible to reduce collision probability and at the same time, as small as possible, to quickly 
select the best relay, before the channel changes again (i.e. within the coherence time of the channel). For example, 
for a mobility of — 3 km/h, the maximum Doppler shift is f m = 2.5 Hz which is equivalent with a minimum 
coherence time on the order of 200 milliseconds. Any relay selection should occur well before that time interval 
with a reasonably small probability of error. From figure 01 we note that selecting c/A w 1/200 will result in a 
collision probability less than 0.6% for policy I. Typical switching times result in c 5/j.s. This gives A « 1ms 
which is two orders of magnitude less than the coherence interval. More sophisticated radios with c as 1/is will 
result in A 200/US, which is three orders of magnitude smaller than the coherence time 6 . 

5 The parenthesized subscripts are due to ordering of the channel gains. 

6 Note that the expected value of the minimum of the set of random variables(timers) is smaller than the average of those random variables. 
So we expect the overhead to be much smaller than the one calculated above 
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A. Calculating Pr(Y 2 < Y\ + c) 

In order to calculate the collision probability from we first need to calculate the joint probability distribution 
of the minimum and second minimum of a collection of M i.i.d 7 random variables, corresponding to the timer 
functions of the M relays. The following theorem provides this joint distribution: 

Theorem 1: The joint probability density function of the minimum and second minimum among M > 2 i.i.d. 
positive random variables T%, T 2 , . . . , Tm, each with probability density function f{t) = dF Jp and cumulative 
distribution function F(t) = Pr(T < t), is given by the following equation: 



fr 1 ,Y 2 (2/1,2/2) 



M (M — 1) f{y x ) f(y 2 ) [1 - F(y 2 )] M - 2 for < Vl < y 2 
elsewhere. 



where Y\ < Y 2 < I3 . . . < Ym are the M ordered random variables T\, T 2 , . . . , Tm- 

Proof: Please refer to appendix U ■ 
Using Theorem 1, we can show the following lemma that gives a closed-form expression for the collision 
probability (equation |9ji: 

Lemma 1: Given M > 2 i.i.d. positive random variables T\, T 2 , . . . , Tm, each with probability density function 
f(x) and cumulative distribution function F(x), and Y\ < Y 2 < I3 . . . < Ym are the M ordered random variables 
T\, T 2 , . . . , Tm, then Pr(Y 2 <Y\ + c), where c > 0, is given by the following equations: 

Pr(Y 2 < Y 1 + c) = l-I c (12) 



/+00 
f(y) [1 - F(y)} M ~ 2 F(y - c) dy (13) 

Proof: 

Please refer to appendix U ■ 
Notice that the statistics of each timer Tj and the statistics of the wireless channel are related according to 
equation Therefore, the above formulation is applicable to any kind of wireless channel distribution. 



B. Results 

In order to exploit theorem [2 and lemma ^ we first need to calculate the probability distribution of Tj for 
i G [1,-M]. From equation © it is easy to see that the cdf F(t) and pdf f(t) of Tj are related to the respective 
distributions of hi according to the following equations: 

F(t)= CDF T% (t) = Pr{T l < t} = 1 - CDF hi (~) (14) 

f(t) = pdf T Xt) = ^-F(t) = A pdf hi (~) (15) 
dt t z t 

After calculating equations (114b . dT5l> . and for a given c calculated from © or ©, and a specific A, we can 
calculate probability of collision using equation (fT2l . 

Before proceeding to special cases, we need to observe that for a given distribution of the wireless channel, 
collision performance depends on the ratio c/A, as can be seen from equation ill Ob . discussed earlier. 

7 The choice of identically distributed timer functions implicitly assumes that the relays are distributed in the same geographical region and 
therefore have similar distances towards source and destination. In that case, randomization among the timers is provided only by fading. 
The cases where the relays are randomly positioned and have in general different distances, is a scenario where randomization is provided 
not only because of fading, but also because of different moments. In such asymmetric cases the collision probability is expected to decrease 
and a concrete example is provided. 
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Fig. 4. Performance in Rayleigh and Ricean fading, for policy I (min) and Policy II (harmonic mean), various values of ratio A/c and 
M = 6 relays, clustered at the same region. Notice that collision probability drops well below 1%. 



1) Rayleigh Fading: Assuming \a S i\, \a,id\, for any i G [1, M], are independent (but not identically distributed) 
Rayleigh random variables, then |a s j| 2 , \aid\ 2 are independent, exponential random variables, with parameters j3±, ftz 
respectively (E[\a si \ 2 } = 1/ftx, E[\a id \ 2 } = l/(3 2 ). 

Using the fact that the minimum of two independent exponential r.v.'s with parameters /?i , /?2 , is again an 
exponential r.v with parameter /3i + fa, we can calculate the distributions for hi under policy I (equation Q. For 
policy II (equation |2j, the distributions of the harmonic mean, have been calculated analytically in [5]. Equations 
(fT4l and ([15} become: 



under policy I: 



F(t) = e~ (l3l+l3 ^ x l l (16) 

m = \(Jh +fr) e-yk+MX/t (17) 
t 2 



under policy II: 



m = A JWW 2 _ x fr+Mm Kl{ WEE.) (18) 



A n a -A (3i+09)/(2t) rPl + A Tjr A V Pi 02 N . n r , A V Pi ft. 



f{t) = __ ftA e-^W)/M [^^^ 1 (^f^ ) + 2^ (^f^ ) ] (19) 
2 i 13 VPl P2 t t 

where Ki(x) is the modified Bessel function of the second kind and order i. 

Equation ( II 2b is calculated for the two policies, for the symmetric case (ft\ = fti = E[\a s i\ 2 ] = E[\aid\ 2 ] = 1) 
of M = 6 relays. Monte-Carlo simulations are also performed under the same assumptions. Results are plotted in 
fig- 13 for various ratios A/c. We can see that Monte-carlo simulations match the results provided by numerical 
calculation of equation ( 112ft with the help of equations (fToT)-(ll9ll. 

Collision probability drops with increasing ratio of A/c as expected. Policy I ("the minimum"), performs 
significantly better than Policy II ("the harmonic mean") and that can be attributed to the fact that the harmonic 
mean smooths the two path SNRs (between source-relay and relay-destination) compared to the minimum function. 
Therefore, the effect of randomization due to fading among the relay timers, becomes less prominent under Policy 
II. The probability can be kept well below 1%, for ratio A/c above 200. 

2) Ricean Fading: It was interesting to examine the performance of opportunistic relay selection, in the case 
of Ricean fading, when there is a dominating communication path between any two communicating points, in 
addition to many reflecting paths and compare it to Rayleigh fading, where there is a large number of equal power, 
independent paths. 
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Assymetry and collision probability 




Keeping the average value of any channel coefficient the same (i?[|a| 2 ] = 1) and assuming a single dominating 
path and a sum of reflecting paths (both terms with equal total power), we plotted the performance of the scheme 
when policy I was used, using Monte-Carlo simulations (fig. |4ji. We can see that in the Ricean case, the collision 
probability slightly increases, since now, the realizations of the wireless paths along different relays are clustered 
around the dominating path and vary less, compared to Rayleigh fading. Policy II performs slightly worse, for the 
same reasons it performed slightly worse in the Rayleigh fading case and the results have been omitted. 

In either cases of wireless fading (Rayleigh or Ricean), the scheme performs reasonably well. 

3 ) Different topologies: For the case of all relays not equidistant to source or destination, we expect the collision 
probability to drop, compared to the equidistant case, since the asymmetry between the two links (from source to 
relay and from relay to destination) or the asymmetry between the expected SNRs among the relays, will increase 
the variance of the timer function, compared to the equidistant case. To demonstrate that, we studied three cases, 
where M = 6 relays are clustered half-way (d/2), closer to transmitter (d/3) or even closer to transmitter (d/W) 
(case 1,2,3 respectively in fig. [5J and d is the distance between source and destination) and one case where the 
relays form an equidistant line network between source and destination (case 4 in fig. |5j). 

Assuming Rayleigh fading, c/A = 1/200 and expected path strength as a non-linear, decreasing function of 
distance (£?[|ay | ] = 1//% oc (l/d^)*), we calculated the collision probability for M = 6 relays, using expressions 
dT5b-dT^b into (11 2t for cases 1, 2, 3 while for case 4 we used Monte-Carlo simulation: in case 1, ft = ft = 1, 
in case 2, ft = (2/3)*, ft = (4/3)* and in case 3, ft = (1/5)*, ft = (9/5)*. For case 4, ft = (2/7)*, ft = 
(12/7)* for the closest terminal to source, ft = (4/7)*, ft = (10/7)* for the second closest terminal to source, 
ft = (6/7)*, ft = (8/7)* for the third closest to source terminal. Due to symmetry, the expected power and 
corresponding f3 factors of the paths, for the third closer to destination, second closer to destination and closest 
terminal to destination, are the same with the ones described before (third closer terminal to source, second closer 
terminal to source and closest to source terminal respectively), with ft and ft interchanged. 

From fig. [5] we can see that the collision probability of asymmetric cases 2, 3 and 4 is strictly smaller compared 
to the symmetric case 1. Policy I performs better than Policy II and collision probability decreases for increasing 
factor v (v = 3, 4 were tested). This observation agrees with intuition that suggests that different moments for the 
path strengths among the relays, increase the randomness of the expiration times among the relays and therefore 
decrease the probability of having two or more timers expire within the same time interval. 

We note that the source can also participate in the process of deciding the best relay. In this special case, 
where the source can receive the CTS message, it could have its own timer start from a value depending upon the 
instantaneous |a s ^| 2 . This will be important if the source is not aware whether there are any relays in the vicinity 
that could potentially cooperate. 
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The proposed method as described above, involving instantaneous SNRs as a starting point for each relay's timer 
and using time (corresponding to an assessment of how good is a particular path within the coherence time of the 
channel) to select space (the best available path towards destination) in a distributed fashion, is novel and has not 
been proposed before, to the best extent of our knowledge. 



IV. Simplifying Cooperative Diversity through Opportunistic Relaying 

We now consider the impact of opportunistic relaying on the cooperative diversity scenario. The main result of this 
section is that opportunistic relaying can be used to simplify a number of cooperative diversity protocols involving 
multiple relays. In particular we focus on the cooperative diversity protocol in [7] which requires the relays to use 
a space-time code while simultaneously transmitting towards the destination. We show that this protocol can be 
simplified considerably by simply selecting the best relay in the second stage. Perhaps surprisingly, this simplified 
protocol achieves the same diversity multiplexing tradeoff achieved in [7]. Furthermore, it does not matter whether 
the relay implements an amplify and forward or a decode and forward protocol in terms of the diversity-multiplexing 
tradeoff. We also note that opportunistic relaying can be used to simplify the non-orthogonal relaying protocols 
proposed in [1]. However the detailed performance analysis is left for future work. 



A. Channel Model 

We consider an i.i.d slow Rayleigh fading channel model following [6]. A half duplex constraint is imposed 
across each relay node, i.e. it cannot transmit and listen simultaneously. We assume that the nodes (transmitter and 
relays) do not exploit the knowledge of the channel at the physical layer. Note that in the process of discovering 
the best relay described in the previous section the nodes do learn about their channel gains to the destination. 
However, we assume that this knowledge of channel gain is limited to the network layer protocol. The knowledge 
of channel gain is not exploited at the physical layer in order to adjust the code rate based on instantaneous channel 
measurements. In practice, the hardware at the physical layer could be quite constrained to allow for this flexibility 
to change the rate on the fly. It could also be that the transmitter is operating at the minimum transmission rate 
allowed by the radio hardware. Throughout this section, we assume that the channel knowledge is not exploited at 
the physical layer at either the transmitter or the relays. 

If the discrete time received signal at the destinationjand the relay node are denoted by Y[n] and Y\ [n] respectively: 
Y[n] = a S( iX[n] + Z[n], n = l,2...— (source transmits destination receives) (20) 

T T 

Y[n] = a r dXi[n] + Z[n], n = —, — + I . . . ,T (best relay transmits destination receives) (21) 

T 

Yi[n] = a sr X[n] + Z\[n\ n = l, 2...— (source transmits best relay receives) (22) 

Here a s d, a r d, a sr are the respective channel gains from the source to destination, best relay to destination 
and source to the best relay respectively. The channel gains between any two pair of nodes are i.i.d M(0, l) 8 . 
The noise Z[n] and Z\[n\ at the destination and relay are both assumed to be i.i.d circularly symmetric complex 
Gaussian jV(0, a 2 ). X[n] and Xi[n] are the transmitted symbols at the transmitter and relay respectively. T denotes 
the duration of time-slots reserved for each message and we assume that the source and the relay each transmit 
orthogonally on half of the time-slots. We impose a power constraint at both the source and the relay: X[n]| 2 ] < P 
and -E[|Xl [n] \ 2 ] < P. For simplicity, we assume that both the source and the relay to have the same power constraint. 
We will define p = P/a 2 to be the effective signal to noise ratio (SNR). This setting can be easily generalized 
when the power at the source and relays is different. 

The following notation is necessary in the subsequent sections of the paper. This notation is along the lines of 
[1] and simplifies the exposition. 

Definition 1: A function f(p) is said to be exponentially equal to b, denoted by f(p) = p b , if 

lim l °f f{p) = b. (23) 

p^co log p 

We can define the relation < in a similar fashion. 

8 The channel gains from the best relay to destination and source to best relay are not JV(0, 1). See Lemma|3|in the Appendix. 
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Definition 2: The exponential order of a random variable X with a non-negative support is given by, 

V = - lim (24) 

p^oo log/9 

The exponential order greatly simplifies the analysis of outage events while deriving the diversity multiplexing 
tradeoff. Some properties of the exponential order are derived in Appendix |H] lemma |2] 

Definition 3: (Diversity-Multiplexing Tradeoff) We use the definition given in [11]. Consider a family of codes 
C p operating at SNR p and having rates R(p) bits per channel use. If P e (R) is the outage probability (see [9]) of 
the channel for rate R, then the multiplexing gain r and diversity order d are defined as 9 

r ^l im ^l ,A log Pe(R) fr) .. 

A-logp d = -hm-j^- (25) 

What remains to be specified is a policy for selecting the best relay. We essentially use the policy 1 (equation 
(0) in the previous section. 

Policy 1: Among all the available relays, denote the relay with the largest value of min{|a sr | 2 , |a r rf|2} as the 
best relay. 

To justify this choice, we note from fig. |3] that the performance of policy I is slightly better than policy II. 
Furthermore, we will see in this section that this choice is optimum in that it enables opportunistic relaying to 
achieve the same diversity multiplexing tradeoff of more complex orthogonal relaying schemes in [7]. We next 
discuss the performance of the amplify and forward and decode and forward protocols. 



B. Digital Relaying - Decode and Forward Protocol 

We will first study the case where the intermediate relays have the ability to decode the received signal, re-encode 
and transmit it to the destination. We will study the protocol proposed in [7] and show that it can be considerably 
simplified through opportunistic relaying. 

The decode and forward algorithm considered in [7] is briefly summarized as follows. In the first half time-slots 
the source transmits and all the relays and receiver nodes listen to this transmission. Thereafter all the relays 
that are successful in decoding the message, re-encode the message using a distributed space-time protocol and 
collaboratively transmit it to the destination. The destination decodes the message at the end of the second time-slot. 
Note that the source does not transmit in the second half time-slots. The main result for the decode and forward 
protocol is given in the following theorem : 

Theorem 2 ([7]): The achievable diversity multiplexing tradeoff for the decode and forward strategy with M 
intermediate relay nodes is given by d(r) = {M + 1)(1 — 2r) for r E (0,0.5). 

The following Theorem shows that opportunistic relaying achieves the same diversity-multiplexing tradeoff if 
the best relay selected according to policy ffl 

Theorem 3: Under opportunistic relaying, the decode and forward protocol with M intermediate relays achieves 
the same diversity multiplexing tradeoff stated in Theorem |2] 

Proof: We follow along the lines of [7]. Let £ denote the event that the relay is successful in decoding 
the message at the end of the first half of transmission and £ denote the event that the relay is not successful in 
decoding the message. Event £ happens when the mutual information between source and best relay drops below 
the code rate. Suppose that we select a code with rate R = r log p and let I(X; Y) denote the mutual information 
between the source and the destination. The probability of outage is given by 

P e = Pt(I(X;Y) < rlogp|£)Pr(£) +Pi(I(X;Y) < rlog p\£)Pv(£) 

= Pr Q log(l + p(\a sd \ 2 + \a rd \ 2 )) < r log p^j Pr(£) + 

Pr Q log(l + p\a sd \ 2 ) < r log p\ Pr(f) 

9 We will assume that the block length of the code is large enough, so that the detection error is arbitrarily small and the main error event 
is due to outage. 
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Pr 


< 


Pr 


< 


P 2r 



< Pr( ^log(l + p(\a sd \ 2 + \a rd \ 2 ))<r\ogp) + 
^log(l + p\a sd \ 2 ) < rlogpj Pr Q 

a sd \ 2 + \a Td \ 2 < p 2r - 1 ) + Pr (\a sd \ 2 < p 2 ^ 1 ) Pr (|a sr | 2 < p 2 ^ 1 
a sd \ 2 < P 27 " 1 ) Pr (M 2 < p 21 '- 1 ) + Pr (|a sd | 2 < p 2 ^ 1 ) Pr (|a sr | z < p lr ~ 

Y p M(2r-\) + p 2r-l p M(2r-l) ^ (Af+l)(2r-l) 

In the last step we have used claim 2 of Lemma |3] in the appendix with m = M. 

■ 

We next study the performance under analog relaying and then mention several remarks. 

C. Analog relaying - Basic Amplify and Forward 

We will now consider the case where the intermediate relays are not able to decode the message, but can only 
scale their received transmission (due to the power constraint) and send it to the destination. 

The basic amplify and forward protocol was studied in [6] for the case of a single relay. The source broadcasts 
the message for first half time-slots. In the second half time-slots the relay simply amplifies the signals it received 
in the first half time-slots. Thus the destination receives two copies of each symbol. One directly from the source 
and the other via the relay. At the end of the transmission, the destination then combines the two copies of each 
symbol through a matched filter. Assuming i.i.d Gaussian codebook, the mutual information between the source 
and the destination can be shown to be [6], ^ 

I(X-Y) = -log (l + p\a sd \ 2 + f{p\a sr \ 2 ,p\a rd \ 2 )) (26) 

/(a, b) = ^— (27) 

The amplify and forward strategy does not generalize in the same manner as the decode and forward strategy 
for the case of multiple relays. We do not gain by having all the relay nodes amplify in the second half of 
the time-slot. This is because at the destination we do not receive a coherent summation of the channel gains 
from the different receivers. If jj is the scaling constant of receiver j, then the received signal will be given by 

y[n] = (y2jLilj a i(ij x [ n ] + z[n]. Since this is simply a linear summation of Gaussian random variables, we do 
not see the diversity gain from the relays. A possible alternative is to have the M relays amplify in a round-robin 
fashion. Each relay transmits only one out of every M symbols in a round robin fashion. This strategy has been 
proposed in [7], but the achievable diversity-multiplexing tradeoff is not analyzed. 

Opportunistic relaying on the other hand provides another possible solution to analog relaying. Only the best 
relay (according to policy is selected for transmission. The following theorem shows that opportunistic relaying 
achieves the same diversity multiplexing tradeoff as that achieved by the (more complicated) decode and forward 
scheme. 

Theorem 4: Opportunistic amplify and forward achieves the same diversity multiplexing tradeoff stated in The- 
orem |2] 

Proof: We begin with the expression for mutual information between the source and destination d26b . An 
outage occurs if this mutual information is less than the code rate r log p. Thus we have that 

P e = Pv(I(X;Y) < rlogp) 

= Pr (log(l + p\a sd \ 2 + /(p|a sr | 2 ,/j|a rd | 2 ) < 2rlogp 

< Pr (\a sd \ 2 < p 2r ~ 1 J{p\a sr \ 2 ,p\a rd \ 2 ) < p 2r ^j 

< Pr (M 2 < p 2r -\min(K r | 2 , |a rd | 2 ) < p 2 ^ 1 + p^^lTp*) 

(b) 

= / 9 2r - 1 p M ( 2r - 1 ) = p (M+l)(2r-l) 
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Fig. 6. The diversity-multiplexing of opportunistic relaying is exactly the same with that of more complex space-time coded protocols. 
Here (a) follows from Lemma |4] and (b) follows from Lemma |3] claim l in appendix |H] and the fact that 

p^V 1 + p 2r -» p 2 ^ 1 as p -» °°- ■ 

D. Discussion 

1) Space-time Coding vs. Relaying Solutions: The (conventional) cooperative diversity setup (e.g. [7]) assumes 
that the cooperating relays use a distributed space-time code to achieve the diversity multiplexing tradeoff in Theorem 
121 Development of practical space-time codes is an active area of research. Recently there has been considerable 
progress towards developing practical codes that achieve the diversity multiplexing tradeoff over MIMO channels. In 
particular it is known that random lattice based codes (LAST) can achieve the entire diversity multiplexing tradeoff 
over MIMO channels [16]. Moreover it is noted in [20] that under certain conditions the analytical criterion such 
as rank and determinant criterion for MIMO links also carry over to cooperative diversity systems 10 . However 
some practical challenges will have to be addressed to use these codes in the distributed antenna setting: (a) The 
codes for MIMO channels assume a fixed number of transmit and receive antennas. In cooperative diversity, the 
number of antennas depends on which relays are successful in decoding and hence is a variable quantity, (b) The 
destination must be informed either explicitly or implicitly which relays are transmitting. 

Opportunistic relaying provides an alternative solution to space time codes for cooperative diversity by using a 
clever relaying protocol. The result of Theorem |4] suggests that there is no loss in diversity multiplexing tradeoff 11 if 
a simple analog relaying based scheme is used in conjunction with opportunistic relaying. Even if the intermediate 
relays are digital, a very simple decode and forward scheme that does eliminates the need for space-time codes can 
be implemented. The relay listens and decodes the message in the first half of the time-slots and repeats the source 
transmission in the second half of the time-slots when the source is not transmitting. The receiver simply does a 
maximal ratio combining of the source and relay transmissions and attempts to decode the message. Theorem [3] 
asserts that once again the combination of this simple physical layer scheme and the smart choice of the relay is 
essentially optimum. 

The diversity-multiplexing tradeoff is plotted in fig. [6j Even though a single terminal with the "best" end-to-end 
channel conditions relays the information, the diversity order in the high SNR regime is on the order of the number 
M + 1 of all participating terminals. Moreover, the tradeoff is exactly the same with that when space-time coding 
across M relays is used. 

2 ) Non-orthogonal Cooperative Diversity Schemes: The focus in this paper was on the multiple relay cooperative 
diversity protocols proposed in [7], since they require that the transmitter and relay operate in orthogonal time-slots 
in addition to the half duplex constraints. The orthogonality assumption was amenable to practical implementation 
[3], since the decoder is extremely simple. More recently, a new class of protocols that relax the assumption that 
the transmitter and relay operate in orthogonal time-slots, (but still assume the half duplex constraint) have been 

'"However it is assumed that the destination knows the channel gain between source and relay for Amplify and Forward, 
"compared to the orthogonal transmission protocols in [7] 
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proposed in [1]. These protocols have a superior performance compared to [7], albeit at the cost of higher complexity 
both at the decoder and network layer. Opportunistic relaying could be naturally used to simplify those protocols 
12 and details of such simplifications and its performance are underway. 

3 ) Impact of Topology: The analysis in for diversity-multiplexing tradeoff was presented assuming that average 
channel gains between each pair of nodes is unity. In other words the impact of topology was not considered. We 
observe that the effect of topology can be included in the analysis using techniques used in [6]. In the high SNR 
regime, we expect fixed multiplicative factors of path loss to contribute little in affecting the diversity-multiplexing 
tradeoff. However topology is certainly important for finite SNR case as observed in [22]. 



V. Conclusion 

We proposed Opportunistic Relaying as a practical scheme for cooperative diversity. The scheme relies on 
distributed path selection considering instantaneous end-to-end wireless channel conditions, facilitates coordination 
among the cooperating terminals with minimum overhead and could simplify the physical layer in communicating 
transceivers by eliminating the requirement of space time codes. 

We presented a method to calculate the performance of the relay selection algorithm, for any kind of wireless 
fading model and showed that successful relay selection could be engineered with reasonable performance. Specific 
examples for Rayleigh and Ricean fading were given. 

We treated Opportunistic Relaying as a distributed virtual antenna array system and analyzed its diversity- 
multiplexing tradeoff, revealing NO performance loss when compared with complex space-time coding protocols 
in the field. 

The approach presented in this work explicitly addresses coordination among the cooperating terminals and has 
similarities with a Medium Access Protocol (MAC) since it directs when a specific node to relay. The algorithm has 
also similarities with a Routing Protocol since it coordinates which node to relay (or not) received information among 
a collection of candidates. Devising wireless systems that dynamically adapt to the wireless channel conditions 
without external means (for example GPS receivers), in a distributed manner, similarly to the ideas presented in 
this work, is an important and fruitful area for future research. 

The simplicity of the technique, allows for immediate implementation in existing radio hardware and its adoption 
could provide for improved flexibility, reliability and efficiency in future 4G wireless systems. 



Appendix I 

Probabilistic Analysis of Successful Path Selection 

Theorem\l\The joint probability density function of the minimum and second minimum among M i.i.d. positive 
random variables X\, X 2 , ■ ■ ■ , Xm, each with probability density function f(x) and cumulative distribution function 
F(x), is given by the following equation: 



fY u Y 2 (yi,y2) 



M (M - 1) f( yi ) f(y 2 ) [1 - F(y 2 )] M - 2 for < y x < y 2 
elsewhere. 



where Y\ < Y 2 < Y$ . . . < Ym are the M ordered random variables X\, X 2 , . . . , Xm- 

Proof: .fa ,Y a (2/1,3/2) dyi dy 2 = Pr(Y 1 G dy 1 , Y 2 G dy 2 ) = 
Pripne X{ in dyi, one Xj in dy 2 (with y 2 > y\ and i ^ j), and all the rest X[s greater than y 2 ) 



= 2 ( A /) Pr( X x G d Vl , X 2 G dy 2 (y 2 > Vl ), X { > y 2 , i G [3,M]) 
= 2 (f )f( yi ) d Vl f(y 2 ) dy 2 [1 - F(y 2 )] M " 2 = 
= M (M — 1) f( yi ) f(y 2 ) [1 - F(y 2 )] M - 2 dy x dy 2 , for < yi < y 2 . 

The third equality is true since there are (^) pairs in a set of M i.i.d. random variables. The factor 2 comes 
from the fact that ordering in each pair matters, hence we have a total number of 2 ( „) cases, with the same 
probability, assuming identically distributed random variables. That concludes the proof. ■ 
Using Theorem 1, we can prove the following lemma: 

12 An Alamouti[15] type code could be used if the relay and source are simultaneously transmitting. 
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Fig. 7. Regions of integration of fY x ,Y 2 (yi, 2/2), for Yi < Y2 needed in Lemma I for calculation of Pr{Yi < Yy + c), c > 0. 



Lemma ^ Given M i.i.d. positive random variables Xl, A" 2 , . . . , Xm. each with probability density function 
f(x) and cumulative distribution function F(x), and Yi < Y 2 < Y3 . . . < Ym the M ordered random variables 
X\, X2, • • • , Xm, then Pr(Y2 < Y\ + c), where c > 0, is given by the following equations: 



Pr(Y 2 <Y x +c) = l-I c 



(28) 



/+00 
/(y) [1 - F(y)] M - 2 F(y - c) dy (29) 

Proof: The joint pdf fy 1; Y 2 (yi > 2/2) integrates to 1 in the region D U -D c , as it can be seen in fig. Q Therefore: 

Pr( Y 2 <Y 1 + c) = J J fofrivuUi) dyl dy2 

I Id ^ Yl ^ m,m ^ dyl dy2 



1 

1 - /, 



Again from fig. I c can easily be calculated: 



+00 

2/2 =c 



J c = m(M-1) / /(y 2 ) [1 - F(y 2 )] JM 



M-2 



3/2 -c 



/(yi) ^yi <iy2 



= M(M-i) / /(y 2 ) [1 - F(y 2 )] M - 2 F(y 2 - c) dy 2 
Jy 2 =c 

The last equation concludes the proof. 

Appendix II 
Diversity-Multiplexing Tradeoff Analysis 

We repeat Definition 1 and Definition 2 in this section for completeness. The relevant lemmas follow. 
Definition 1: A function f(p) is said to be exponentially equal to b, denoted by f(p) = p b , if 

lim = ft. 

p-^co log p 



(30) 



(31) 



We can define the relation < in a similar fashion. Definition 2: The exponential order of a random variable X 
with a non-negative support is given by, 

log X 

(32) 



V 



lim 

p^QO log p 



Lemma 2: Suppose Xi, X 2 , • • • , X m are m i.i.d exponential random variables with parameter A (mean 1/A), and 
X = max{Xi, X 2 , ■ . . X m }. If V is the exponential order of X then the density function of V is given by 



fv{v) 



v>0 
v < 



(33) 
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and 

Pr(X < p~ v ) = p- mv (34) 

Proof: Define, 

V =-**±. 

Thus V p is obtained from definition |2j without the limit of p — ► oo. 

Pr(Vp >v) = Pr(X < p- v ) 

= Pr(Xx < X 2 < . . . X m < p- v ) 

m 

= l[PT(Xi<p- V ) 

i=l 

= (l-exp(-A^)) m 



Note that Pr(V r p > v) m p mv . Differentiating with respect to v and then taking the limit p — ► cxd, we recover 

EJ. ■ 

From the above it can be seen that for the simple case of a single exponential random variable (m = 1), 
Pi(X < p- v ) = Pr(Vp >v) = p- v . 

Lemma 3: For relays, j = 1, 2, . . . , m, let a s j and ajd denote the channel gains from source to relay j and relay 
j to destination. Suppose that a sr and aw denote the channel gain of the source to the best relay and the best relay 
to the destination, where the relay is chosen according to rule ^ i.e. 

min(|a sr | 2 , |a rrf | 2 ) = max{min(|a s i| 2 , \a ld \ 2 ), . . . ,min(|a sm | 2 , |a m d| 2 )} 

Then, 

1) min(|a sr | 2 , |ow| 2 ) has an exponential order given by d33l . 

p -mv y > o 
1 otherwise 



Pr(K r | 2 <p~ v ) =Pr(|a rd | 2 < p~ v ) < 



Proof: Let us denote = min(|a s j| 2 , la^l 2 )- Since each of the X^ are exponential random variables with 
parameter 2, claim 1 follows from Lemma|2] Also since \a s d\ 2 and \a rd \ 2 cannot be less than min(|a s rf| 2 , |a r ^| 2 ) 
claim 2 follows immediately from claim 1. ■ 
Lemma 4: With /(•,•) defined by relation d27l . we have that 

Pr (f(pa,pb) < p 2r ) < Pr (mm(a,b) < p 2r ~ x + p r - 1 ^/TTr 2 ' 
Proof: Without loss in generality, assume that a > b. 

f(pa,pb) = p ~T~TT 
a + b+- 

= pb 



a + h+ \, 



(a) 

> pb 



26+ i 



p / 

Here (a) follows since is an increasing function in a, for K > and a>b. 
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Now we have that 



Pr(f(pa,pb)<p 2r ) < Pr^-J^Kp^-i 



(a) 

Where (a) follows since b > so that Pr(6 < 0) = 0. 



Pr(6 2 <2p 2r - 1 b + p 2r - 2 ) 
Pi((b-p 2r - 1 )2<p* r - 2 + p 2r - 2 ) 

Pr (b < p 2r - x + p r - x ^fl~Vp^) 
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